Decoherence : An Irreversible Process 



Roland OmnesQ 

Laboratoire de Physique Theoriqu^ 
Universite de Paris XI, Bdtiment 210, 91405 Orsay Cedex, France 



Abstract 

A wide-ranging theory of decoherence is derived from the quantum theory 
of irreversible processes, with specific results having for their main limitation 
the assumption of an exact pointer basis. 
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Decoherence has become widely recognized as an essential step for understanding 
and interpreting quantum mechanics. It has been mainly investigated on solvable 
models, which give most of what we know about the effect ; but though some 
of these models are considered as realistic in a definite situation (for instance in 
quantum optics [jl|), they remain in most cases rather far from reality. Now that 
the importance of decoherence is acknowledged, it certainly requires a wider theory, 
and forthcoming investigations will also need more precise results (particularly in 
the research on quantum computers, where decoherence is expected to be the main 
obstacle and must be controlled quantitatively). 

A simple starting point for such a theory consists in noting that decoherence is 
an irreversible process, so that one can apply the rather general theories for such 
processes. It is generally agreed that the best theories of that sort rely on the so- 
called "projection method" [0-0]. There have already been some attempts to apply 
it to decoherence 0, but the results in the present Letter improve on them on at 
least three significant points : 

1. One can always define an average part in the coupling of the system with 
the environment, and when that part is removed the remaining coupling con- 
sists only of fluctuations ; one can therefore extend considerably the range of 
perturbation theory. 

2. A rather general master equation for decoherence is obtained (Eq. (|12]) below). 

3. Precise quantitative equations are obtained for decoherence, at least when an 
exact "pointer basis" exists |^. 



As far as the overall framework is concerned, one assumes as usual that the 
system under consideration can be split into a "collective" subsystem with hamil- 
tonian and an environment having a very large number of degrees of freedom 
with hamiltonian He- The two systems are coupled, and the full hamiltonian is 
H = He® le + Ic® He + Hi, the coupling Hi allowing energy exchanges and other 
mutual influences (including decoherence) between the collective system and the 
environment. The full density operator p evolves according to the basic equation 
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P = -i[H,p] 



(1) 



The main features of the projection method for a typical irreversible process 
are as follows. A (countable or not) set of independent "relevant observables" 
(including the identity /) is selected. Their "exact" average values, resulting from 
the exact density operator p{t), are denoted by a*(t). A time-dependent test density 
operator po{t) is then introduced with the assumptions (i) that it gives the exact 
average values {a*(t)} for the relevant operators {^4*} (so that one might know them 
if po is known), and (ii) that its information content is minimal. It must be of the 
form 

Po = exp (-A, A') , (2) 

with Lagrange parameters Aj ; summation over repeated indices is assumed as usual, 
and Po is normalized since the identity / belongs to the set of relevant observables. 
Auxiliary "densities" (or more properly trace-class operators) are then defined by 
Si = dpo/da\ They satisfy the orthogonality properties, 

Tr{siA^)=6l , (3) 

which amount essentially to da^/da^ = 6j. 

The theory makes use of "superoperators" , acting linearly on a trace-class oper- 
ator to yield a similar operator. For instance, Eq. (|I]) can be written conventionally 
as p = £p, where C is the so-called Liouville superoperator. Another superoperator 
is defined in the projection approach by 

V = s,^A , (4) 

which means that V acts on a trace-class operator p to give Vfi = Tr{A^ p)si. Eq. 
d^) imphes the projection property = V. A relevant density operator is defined 
as pi = Vp. It also yields the exact quantities {a*} as average values of the relevant 
operators {A*}, in view of Eq. (H). Denoting by J the identity superoperator, one 
also introduces the superoperator Q = J — V which satisfies the same projection 
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property = Q. Denoting Qp by p2 (so that p = pi + P2), and applying V and Q 
to both sides of Eq. (|T|) one obtains evolution equations for pi and p2 : 

pi = VCVpi + VVpi + VCQp2 + VQp2 , (5) 
P2 = QCQp2 -VQp2 + QCVpi - VVpi . (6) 

When applying the projection method to decoherence, it will be convenient to 
introduce a commuting set of collective observables X whose eigenvalues x are either 
discrete or continuous. The "relevant observables" {A^} are chosen to consist of 
the identity /, the environment hamiltonian Ic ® He and the collective observables 
(|x > ±\x' >) ■ (< x|± < x'l) ® Je and (|x > ±i\x' >) ■ (< x| =F ^ < x') ® Jg for every 
pair (x, x') of eigenvalues of X. One can use more simply the set of non-hermitian 
operators A^^' = \x >< x'\ ® le, which clearly provide a basis for the collective 
observable. The test density operator (|^) takes then the simple form 

po = pc® Pe , (7) 

where pc turns out to be the familiar reduced density operator for the collective 
subsystem and pe is, at least formally, a normalized density operator for the envi- 
ronment as if it were in thermal equilibrium : 

Pc = tr p , pe = exp {-a - l3He) . (8) 

Throughout I denote a partial trace over the environment by tr and a full trace by 
Tr. The time-depending parameters a and (3 are chosen so that pe is normalized 
and the "exact" average value E for the environment energy is obtained from it. It 
should be stressed that this expression of pe does not mean that the environment is 
in thermal equilibrium ; it means only that one does not need to know more than 
the average environment energy in order to obtain collective quantities, including 
the reduced density operator. 

The first part of the calculation consists in obtaining algebraically obtaining the 
auxiliary densities Sj, which are given here for convenience (using the notation s(A') 
in place of Si) 
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S {\X >< X'\ le) = \x' >< X'\ (g) Pe , 

s{h0H,)=h®{p,{H,-E))A-^ ; s{I) = -E h ® {p, {H, - E)) A'^ 

with = tr{Hlpe) — E'^, so that acting on a trace-class operator p, the projection 
superoperator V gives 

Vp = [tr p)®pe + Ic® {pe {He - E)) A"" {Tt {H^ p)~ETr p} . (9) 

(In particular, pi = po). 

Then comes an important trick. It is very convenient to introduce an average 
collective coupling 

AHc = tr {Hi ■ Je ® Pe) , 

which is a collective operator representing a collective effect of the environment 
(for instance the action of pressure in the case of a gaseous environment). It is 
generally important, though equal to zero in a few special cases (matter-radiation 
coupling, nuclear magnetic resonance, and some oscillator models). The remaining 
part of the coupling H[ = Hi — AHe ® /g, consists presumably in most cases of 
small fluctuations, which can be considered as perturbations of the hamiltonian 
Hq = {He + AHc) ® le + Ic ® He- To sccoud ordcr in H'l, Eqs. (|§-|^) become 
explicitly 

Pc = -ilHc, Pc] - i tr[Hi,p2] , (10) 

P2 = -i[Ho, P2] - i[Hi,po] + i tr[Ho, P2] pe ■ (H) 

One can solve Eq. (^) for p2 as a function of po, using first-order perturbation theory 
and assuming for convenience that the environmement is in thermal equilibrium at 
an initial time t = so that p2(0) = 0. Inserting the resulting expression for p2 into 
Eq. ([ToD , one obtains the rather simple and fundamental "master equation" 
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t[He + AHe, Pc] - dt' tr { [h'i, U[H'i,pQ{t')]U-^] } , (12) 
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where U = exp{—iHQ{t — t')). 



Though aheady known when the full coupling Hi is weak, Eq. (p!2|) has a much 
wider range of valid since it holds for a fluctuating H[ The last step in the cal- 
culation consists in writing down explicitly the master equation, explicit expressions 
being obtained when the basis |x > is an exact pointer basis or, more precisely, when 
H[ and X commute so that H[ is diagonal in the \x > basis and behaves like an 
operator V{X) in the environment Hilbert space (more explicitly, introducing eigen- 
vectors \k > of He with eigenvalues Ek, one has < x, k\H[\x', n >= 6{x — x')Vkn{x)). 
The existence of a pointer basis will be assumed from here on. 

One can then define a microscopic distance (abbreviated by fiD) between two 
points X and distance |x — x'| where the (quantum) first term in the right- 

hand side of Eq. (|12|) dominates the value of Pc{x, x'). A small macroscopic distance 
(SMD) will be one for which the second (decoherence) term in the right-hand side 
dominates, although x — x' is still macroscopically small. I will assume furthermore 
that V{x) — V{x') depends linearly locally on x — x' with a "slope" V over small 
distances {fiD and SMD), which are the only distances of interest in applications. 
The second (decoherence) term in the right-hand side of the master equation (|I2|) 
becomes then, for fiD^s and SMD^s, 

— J dt' J dxi dx[ K{x, x'; Xi, x[;t — t')pc{xi, x[,t') (13) 
where the kernel is 



K = ^ (x — x') [{xi — x'l) cosh Puknf^ + {xi + x[ — X — x') sinh jSukn/'A 

n,k,N,M 

X < x\N >< N\xi >< x[\M >< M\x' > exp{iilMNT) 
x\V^'^exp{-iuJknT)pnk ■ (14) 

By >, \M >, one denotes eigenvectors of H^ + AHc with eigenvalues E^, Em, 
and Qmn = Em - En, similarly Ukn = Ek-En and pnk = exp[-a -(3{Ek + En)/2), 
and T = t — t' . 

This expression has many interesting consequences : 
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1. The results of previous models can be recovered, often resulting in a simpler 
kernel ; for instance, models with an environment consisting of a collection 
of two-states systems 0] (one system for each u), or harmonic oscillators 0. 
Decoherence by collisions with an environment of molecules or photons 
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best obtained by using plane-wave states for n and outgoing scattering states 
for k. 

2. Both terms involving cosh/5a;/2 and sinh/5u;/2 in Eq. ( |14|) are significant at 
low temperature and fiD^s. This case will be presumably important for future 
technology when decoherence and quantum coherence compete. 

3. The interpretation of quantum measurements, with suppression of macroscopic 
superpositions (as in the Schrodinger cat problem), is mainly concerned with 
macroscopic values of (x — x'), or SMD^s. One can then put Xi = x and x'l = 
Xi, as can be shown easily when = P"^ /2m by means of Fourier transforms. 
The case He = P'^ /2m + W{x) requires a more elaborate justification using 
coherent collective states or microlocal analysis, but one always obtains for 
SMD's 



K ~ 5{x — xi)5{x' — x\){x — x'Y cosh.{l3uJkn/2) 

n,k 

x|Vfc'„|^exp[-ZC^fcn(t-t')]Pnfc (15) 

4. At high enough temperature and when retardation m t — t' is neglected, the 
decoherence term at SMD's becomes simply —fi{x — x')^, with a decoherence 
coefficient 

- - II \V'L\K^kn - iOy^Pnk - TT ^ I PPnfc5(Wfcn) (16) 
n,k n,k 

5. If P is the momentum canonically conjugate to X, one finds easily that the 
term in cosh/5u;/2 in Eq. (p^lf ) does not contribute to < dP/dt >, i.e. to 
damping. The term in sinh (3uj / 2 gives on the other hand gives a damping 

< dP/dt > = forces - /* D(t - t') < Pit') > dt' , (17) 

Jo 
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with 



D{t) = 2m-i 



'kn 



This simple result results from several steps : The action of P = —id/dx on 
the kernel K (in Eq. (|I^) ) removes the factor (x — x') since other terms in 
the derivative containing this factor have zero average. The integral over x of 
< x\N >< M\x > gives simply 6nm- A term such as x < a;|A^ >< M\x > 
exp(zfir) is expressed as < M\U^XUc\N >, with Uc = exp{—iHcT). One can 
then use the equation 



which is exact for He = P^/2m and valid up to higher orders in h in the 
presence of a collective potential. The factor r in Eq. ([ISD is finally removed by 
integrating by parts on r. Other more direct methods also exist for evaluating 
damping, and they agree with the result. 

When retardation effects are again neglected, the damping term in Eq. 
(0) becomes —7 < P > and, at high enough temperature, one recovers the 
well-known relation fx = mkTh~'^, where Planck's constant has been reintro- 
duced. One might also derive the familiar expression —7/2(0; — x'){d/dx — 
d/dx')pc{x,x') for damping with similar approximations. 

To conclude, the main limitation of the final results following Eq. (|12D, is the 
assumption of an exact pointer basis. Such bases are known to exist for a strictly me- 
chanical collective system (the X's being the position coordinates of coarse-grained 
pieces of matter in the macroscopic system), or in the case of SQUID loops |^. More 
generally, approximate diagonalization occurs also presumably in a basis of coherent 
states ||ll|] although the validity of the present analysis remains uncertain in that 
case, as well as its relation with classical behavior il^l . 

Finally, the increase of entropy resulting from irreversibility is mainly found in 
the increase of —Tr{pc Log pc) when pc becomes approximately diagonal. 
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[X, U] = UPr/m 



(18) 
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